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The ATLAS experiment has measured the production cross-section of events with two isolated 
photons in the final state, in proton-proton collisions at ^/s = 7 TeV. The full data set acquired in 
2010 is used, corresponding to an integrated luminosity of 37 pb~^. The background, consisting of 
hadronic jets and isolated electrons, is estimated with fully data-driven techniques and subtracted. 
The differential cross-sections, as functions of the di-photon mass {m-y-y), total transverse momentum 
{Pt,-ii) s-'^d azimuthal separation (At/)^^,), are presented and compared to the predictions of next- 
to-leading-order QCD. 



I. INTRODUCTION 

The production of di-photon final states in proton- 
proton collisions may occur through quark- antiquark t- 
channel annihilation, qq ^ 77, or via gluon-gluon inter- 
actions, gg — ?> 77, mediated by a quark box diagram. 
Despite the higher order of the latter, the two contribu- 
tions are comparable, due to the large gluon flux at the 
LHC. Photon-parton production with photon radiation 
also contributes in processes such as qq, gg — >■ 577 and 
qg — >■ (777. During the parton fragmentation process, 
more photons may also be produced. In this analysis, all 
such photons are considered as signal if they are isolated 
from other activity in the event. Photons produced after 
the hadronization by neutral hadron decays, or coming 
from radiative decays of other particles, are considered 
as part of the background. 

The measurement of the di-photon production cross- 
section at the LHC is of great interest as a probe of 
QCD, especially in some particular kinematic regions. 
For instance, the distribution of the azimuthal separa- 
tion, A(J3-yy , is sensitive to the fragmentation model, espe- 
cially when both photons originate from fragmentation. 
On the other hand, for balanced back-to-back di-photons 
{A(f>jj ~ TT and small total transverse momentum, PT.77) 
the production is sensitive to soft gluon emission, which 
is not accurately described by fixed-order perturbation 
theory. 

Di-photon production is also an irreducible background 
for some new physics processes, such as the Higgs decay 
into photon pairs [ij: in this case, the spectrum of the 
invariant mass, rrij-y, of the pair is analysed, searching for 
a resonance. Moreover, di-photon production is a charac- 
teristic signature of some exotic models beyond the Stan- 
dard Model. For instance. Universal Extra-Dimensions 
(UED) predict non-resonant di-photon production asso- 
ciated with significant missing transverse energy [2|, |3|. 
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Other extra-dimension models, such as Randall-Sundrum 
[J], predict the production of gravitons, which would de- 
cay into photon pairs with a narrow width. 

Recent cross-section measurements of di-photon pro- 
duction at hadron colliders have been performed by the 
D0 Q and CDF jS] collaborations, at the Tevatron 
proton-antiproton collider with a centre-of-mass energy 
yi = L96 TeV. 

In this document, di-photon production is studied in 
proton-proton collisions at the LHC, with a centre-of- 
mass energy y/s = 7 TeV. After a short description of 
the ATLAS detector (Section [H]), the analysed collision 
data and the event selection are detailed in Section IIIIl 
while the supporting simulation samples are listed in Sec- 
tion |lVl The isolation properties of the signal and of the 
hadronic background are studied in Section |Vl The eval- 
uation of the di-photon signal yield is obtained by sub- 
tracting the backgrounds from hadronic jets and from 
isolated electrons, estimated with data-driven methods 
as explained in Section IVll Section |VlTl describes how the 
event selection efficiency is evaluated and how the final 
yield is obtained. Finally, in Section IVIIIl the differen- 
tial cross-section of di-photon production is presented as 
a function of m.^^, pt,77 and A(/)^^. 



II. THE ATLAS DETECTOR 

The ATLAS detector [7| is a multipurpose particle 
physics apparatus with a forward-backward symmetric 
cylindrical geometry and near Att coverage in solid angle. 
ATLAS uses a right-handed coordinate system with its 
origin at the nominal interaction point (IP) in the cen- 
tre of the detector and the z-axis along the beam pipe. 
The X-axis points from the IP to the centre of the LHC 
ring, and the y axis points upward. Cylindrical coor- 
dinates (r, 0) are used in the transverse plane, (j) being 
the azimuthal angle around the beam pipe. The pseu- 
dorapidity is defined in terms of the polar angle 9 as 
rj = — lntan(0/2). The transverse momentum is defined 



as pt = psinO = p/ cosh rj, and a similar definition holds 
for the transverse energy Ej^. 

The inner tracking detector (ID) covers the pseudo- 
rapidity range \r]\ < 2.5, and consists of a silicon pixel 
detector, a silicon microstrip detector, and a transition 
radiation tracker in the range \r]\ < 2.0. The ID is sur- 
rounded by a superconducting solenoid providing a 2 T 
magnetic field. The inner detector allows an accurate 
reconstruction of tracks from the primary proton-proton 
collision region, and also identifies tracks from secondary 
vertices, permitting the efficient reconstruction of pho- 
ton conversions in the inner detector up to a radius of 
« 80 cm. 

The electromagnetic calorimeter (ECAL) is a lead- 
liquid argon (LAr) sampling calorimeter with an accor- 
dion geometry. It is divided into a barrel section, cov- 
ering the pseudorapidity region \ri\ < 1.475, and two 
end-cap sections, covering the pseudorapidity regions 
1.375 < |?7| < 3.2. It consists of three longitudinal layers. 
The first one, in the ranges \ri\ < 1.4 and 1.5 < \ri\ < 2.4, 
is segmented into high granularity "strips" in the 77 di- 
rection, sufficient to provide an event-by-event discrimi- 
nation between single photon showers and two overlap- 
ping showers coming from a 7r° decay. The second layer 
of the electromagnetic calorimeter, which collects most 
of the energy deposited in the calorimeter by the photon 
shower, has a thickness of about 17 radiation lengths and 
a granularity of 0.025 x 0.025 in 77 x </> (corresponding to 
one cell). A third layer is used to correct leakage beyond 
the ECAL for high energy showers. In front of the ac- 
cordion calorimeter a thin presampler layer, covering the 
pseudorapidity interval jr/j < 1.8, is used to correct for 
energy loss before the calorimeter. 

The hadronic calorimeter (HCAL), surrounding the 
ECAL, consists of an iron-scintillator tile calorimeter 
in the range {rjl < 1.7, and two copper-LAr calorime- 
ters spanning 1.5 < |ry| < 3.2. The acceptance is ex- 
tended by two tungsten-LAr forward calorimeters up to 
|?7| < 4.9. The muon spectrometer, located beyond the 
calorimeters, consists of three large air-core supercon- 
ducting toroid systems, precision tracking chambers pro- 
viding accurate muon tracking over \r]\ < 2.7, and fast 
detectors for triggering over \r]\ < 2.4. 

A three-level trigger system is used to select events 
containing two photon candidates. The first level trigger 
(level-1) is hardware based: using a coarser cell granu- 
larity (0.1 X 0.1 in ?7 X 0), it searches for electromagnetic 
deposits with a transverse energy above a programmable 
threshold. The second and third level triggers (collec- 
tively referred to as the "high-level" trigger) are imple- 
mented in software and exploit the full granularity and 
energy calibration of the calorimeter. 



III. COLLISION DATA AND SELECTIONS 

The analysed data set consists of proton-proton colli- 
sion data at -y/s = 7 TeV collected in 2010, corresponding 



to an integrated luminosity of 37.2 ± 1.3 pb~^ Q. The 
events are considered only when the beam condition is 
stable and the trigger system, the tracking devices and 
the calorimeters are operational. 



A. Photon reconstruction 

A photon is defined starting from a cluster in the 
ECAL. If there are no tracks pointing to the cluster, the 
object is classified as an unconverted photon. In case of 
converted photons, one or two tracks may be associated 
to the cluster, therefore creating an ambiguity in the clas- 
sification with respect to electrons. This is addressed as 
described in Ref [3] . 

A fiducial acceptance is required in pseudorapidity, 
\ri''\ < 2.37, with the exclusion of the barrel/endcap tran- 
sition 1.37 < \ri''\ < 1.52. This corresponds to the regions 
where the ECAL "strips" granularity is more effective for 
photon identification and jet background rejection [9|. 
Moreover, photons reconstructed near to regions affected 
by read-out or high- voltage failures are not considered. 

In the considered acceptance range, the uncertainty 
on the photon energy scale is estimated to be ^ ±1%. 
The energy resolution is parametrized as ue/E ~ 
a/ ^JE [GeV] ® c, where the sampling term a varies be- 
tween 10% and 20% depending on rf , and the constant 
term c is estimated to be 1.1% in the barrel and 1.8% 
in the endcap. Such a performance has been measured 
in Z — >■ e'^e~ events observed in proton-proton collision 
data in 2010. 



B. Photon selection 

The photon sample suffers from a major background 
due to hadronic jets, which generally produce calori- 
metric deposits broader and less isolated than electro- 
magnetic showers, with sizable energy leaking to the 
HCAL. Most of the background is reduced by apply- 
ing requirements (referred to as the LOOSE selection, L) 
on the energy fraction measured in the HCAL, and on 
the shower width measured by the second layer of the 
ECAL. The remaining background is mostly due to pho- 
ton pairs from neutral hadron decays (mainly tt*^) with 
a small opening angle, and reconstructed as single pho- 
tons. This background is further reduced by applying a 
more stringent selection on the shower width in the sec- 
ond ECAL layer, together with additional requirements 
on the shower shape measured by the first ECAL layer: 
a narrow shower width and the absence of a second sig- 
nificant maximum in the energy deposited in contiguous 
strips. The combination of all these requirements is re- 
ferred to as the tight selection (T). Since converted 
photons tend to have broader shower shapes than uncon- 
verted ones, the cuts of the TIGHT selection are tuned 
differently for the two photon categories. More details 
on these selection criteria are given in Ref jlOj . 



To reduce the jet background further, an isolation 
requirement is apphed: the isolation transverse energy 
E^°, measured by the calorimeters in a cone of angular 
radius R = ^{rj — rfY + (</*" 4''^Y < 0-4, is required to 
satisfy El^° < 3 GeV (isolated photon, I). The calcula- 
tion of i?^° is performed summing over ECAL and HCAL 
cells surrounding the photon candidate, after removing a 
central core that contains most of the photon energy. An 
out-of-core energy correction [10| is applied, to make E^° 
essentially independent of E^, and an ambient energy 
correction, based on the measurement of soft jets 11, ij] 
is applied on an event- by- event basis, to remove the con- 
tribution from the underlying event and from additional 
proton-proton interactions ( "in-time pile-up" ) . 



C. Event selection 

The di-photon candidate events are selected according 
to the following steps: 

• The events are selected by a di-photon trigger, 
in which both photon candidates must satisfy 
the trigger selection and have a transverse energy 
E^ > 15 GeV. To select genuine collisions, at least 
one primary vertex with three or more tracks must 
be reconstructed. 

• The event must contain at least two photon candi- 
dates, with E^ > 16 GeV, in the acceptance defined 
in Section IIII Al and passing the LOOSE selection. 
If more than two such photons exist, the two with 
highest E^ are chosen. 

• To avoid a too large overlap between the two 
isolation cones, an angular separation AR^j = 



VirPi 



:'( 



vl? 



b\Y > 0.4 is required. 



• Both photons must satisfy the tight selection (TT 
sample). 

• Both photons must satisfy the isolation require- 
ment E}^° < 3 GeV (TITI sample). 

In the analysed data set, there are 63673 events where 
both photons satisfy the LOOSE selection and the AR^^ 
separation requirement. Among these, 5365 events be- 
long to the TT sample, and 2022 to the TITI sample. 



IV. SIMULATED EVENTS 

The characteristics of the signal and background events 
are investigated with Monte Carlo samples, generated us- 
ing Pythia 6.4.21 [13]. The simulated samples are gener- 
ated with pile-up conditions similar to those under which 
most of the data were taken. Particle interactions with 
the detector materials are modelled with Geant4 [ij| 
and the detector response is simulated. The events are 
reconstructed with the same algorithms used for collision 



data. More details on the event generation and simula- 
tion infrastructure are provided in Ref |15| . 

The di-photon signal is generated with Pythia, 
where photons from both hard scattering and quark 
bremsstrahlung are modelled. To study systematic ef- 
fects due to the generator model, an alternative di-photon 
sample has been produced with Sherpa fl6]. 

The background processes are generated with the main 
physical processes that produce (at least) two sizable 
calorimetric deposits: these include di-jet and photon- 
jet final states, but minor contributions, e.g. from W, Z 
bosons, are also present. Such a Monte Carlo sample, 
referred to as "di-jet-like" , provides a realistic mixture 
of the main final states expected to contribute to the 
selected data sample. Moreover, dedicated samples of 
W ^ ev and Z — >■ e'^e" simulated events are used for 
the electron/photon comparison in isolation and back- 
ground studies. 



V. PROPERTIES OF THE ISOLATION 
TRANSVERSE ENERGY 

The isolation transverse energy, _E^° , is a powerful dis- 
criminating variable to estimate the jet background con- 
tamination in the sample of photon candidates. The ad- 
vantage of using this quantity is that its distribution can 
be extracted directly from the observed collision data, 
both for the signal and the background, without relying 
on simulations. 

Section IV Al describes a method to extract the distri- 
bution of i?^° for background and signal, from observed 
photon candidates. An independent method to extract 
the signal E!^° distribution, based on observed electrons, 
is described in Section IV Bl Finally, the correlation be- 
tween isolation energies in events with two photon can- 
didates is discussed in Section IV CI 



Background and signal isolation from photon 
candidates 



For the background study, a control sample is defined 
by reconstructed photons that fail the tight selection 
but pass a looser one, where some cuts are released on 
the shower shapes measured by the ECAL "strips" . Such 
photons are referred to as non-TiGHT. A study carried 
out on the "di-jet-like" Monte Carlo sample shows that 
the E}^° distribution in the non-TiGHT sample reproduces 
that of the background, as shown in Figure [IJa). 

The TIGHT photon sample contains a mixture of sig- 
nal and background. However, a comparison between 
the shapes of the E!^° distributions from tight and non- 
TIGHT samples (FigureJUb)) shows that for i?^° > 7 GeV 
there is essentially no signal in the tight sample. There- 
fore, the background contamination in the tight sample 
can be subtracted by using the non-TiGHT sample, nor- 
malized such that the integrals of the two distributions 
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FIG. 1. Extraction of the isolation energy {E^°) distribu- 
tions, for signal and background. The plots are made with 
a "di-jet-like" Monte Carlo sample: the "signal" and "back- 
ground" classifications are based on the Monte Carlo infor- 
mation, (a) Normalized E^° distribution for the background 
and for the non-TlGHT sample, (b) E^° distribution, for the 
TIGHT and the non-TlGHT samples: the latter is scaled as ex- 
plained in the text, (c) Normalized E^° distribution for the 
signal and for the tight sample, after subtracting the scaled 
non-TIGHT sample. In (a,c) the vertical line shows the isola- 
tion cut £;^=° < 3 GeV. 



are equal for Eif° > 7 GeV. The E}§° distribution of the 
signal alone is thus extracted. Figure ^c) shows the re- 
sult, for photons in the "di-jet-like" Monte Carlo sample. 

In collision data, events with two photon candidates 
are used to build the tight and non-TiGHT samples, for 
the leading and subleading candidate separately. The 
points in Figure [2] display the distribution of i?^° for 
the leading and sub-leading photons. In each of the two 
distributions, one bin has higher content, reflecting op- 
posing fluctuations in the subtracted input distributions 
in those bins. The effect on the di-photon cross-section 
measurement is negligible. 



The main source of systematic error comes from the 
definition of the non-TiGHT control sample. There are 
three sets of strips cuts that could be released: the first 
set concerns the shower width in the core, the second 
tests for the presence of two maxima in the cluster, and 
the third is a cut on the full shower width in the strips. 
The choice adopted is to release only the first two sets 
of cuts, as the best compromise between maximizing the 
statistics in the control sample, while keeping the back- 
ground E}^° distribution fairly unbiased. To test the ef- 
fect of this choice, the sets of released cuts have been 
changed, either by releasing only the cuts on the shower 
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FIG. 2. Data-driven signal isolation distributions for the lead- 
ing (top) and sub-leading (bottom) photons obtained using 
the photon candidates (solid circles) or extrapolated from 
electrons (continuous lines). 



core width in the strips, or by releasing all the strips 
cuts. A minor effect is also due to the choice of the re- 
gion i?^° > 7 GeV, to normalize the non-TIGHT control 
sample: the cut has therefore been moved to 6 and 8 GeV. 
More studies with the "di-jet-like" Monte Carlo sam- 
ple have been performed, to test the robustness of the 
E}^° extraction against model-dependent effects such as: 
(i) signal leakage into the non- tight sample; (ii) corre- 
lations between E!^° and strips cuts; (Hi) different sig- 
nal composition, i.e. fraction of photons produced by 
the hard scattering or by the fragmentation process; (iv) 
different background composition, i.e. fraction of pho- 
ton pairs from tt" decays. In all cases, the overall sys- 
tematic error, computed as described above, covers the 
differences between the true and data-driven results as 
evaluated from these Monte Carlo tests. 



B. Signal isolation from electron extrapolation 

An independent method of extracting the Ei^° distri- 
bution for the signal photons is provided by the "electron 



extrapolation" . In contrast to photons, it is easy to select 
a pure electron sample from data, from W -^ e^v and 
Z -^ e+e~ events [l^l- The main differences between the 
electron and photon £'^° distributions are: (i) the elec- 
tron E}^° in the bulk of the distribution is slightly larger, 
because of bremsstrahlung in the material upstream of 
the calorimeter; {ii) the photon ElfP distribution exhibits 
a larger tail because of the contribution of the photons 
from fragmentation, especially for the sub-leading pho- 
ton. Such differences are quantified with W^ — >■ e^u, 
Z — > e+e^ and 77 Monte Carlo samples by fitting the 
E}f distributions with Crystal Ball functions [la] and 
comparing the parameters. Then, the electron/photon 
differences are propagated to the selected electrons from 
collision data. The result is shown by the continuous 
lines in Figure [H agreeing well with the ElfP distribu- 
tions obtained from the non-TiGHT sample subtraction 
(circles). 



C. Signal and background isolation in events with 
tw^o photon candidates 

In events with two photon candidates, possible cor- 
relations between the two isolation energies have been 
investigated by studying the signal and background El^° 
distributions of a candidate ( "probe" ) under different iso- 
lation conditions of the other candidate ("tag"). The 
signal ElfP shows negligible dependence on the tag con- 
ditions. In contrast, the background ElfP exhibits a clear 
positive correlation with the isolation transverse energy 
of the tag: if the tag passes (or fails) the isolation require- 
ment, the probe background candidate is more (or less) 
isolated. This effect is visible especially in di-jet final 
states, which can be directly studied in collision data by 
requiring both photon candidates to be non-TiGHT, and 
is taken into account in the jet background estimation 
(Section EEl). 

This correlation is also visible in the "di-jet-like" 
Monte Carlo sample. 



VI. BACKGROUND SUBTRACTION AND 
SIGNAL YIELD DETERMINATION 

The main background to selected photon candidates 
consists of hadronic jets. This is reduced by the pho- 
ton TIGHT selection described in Section fill Bl However 
a significant component is still present and must be sub- 
tracted. The techniques to achieve this are described in 

Section EEl 

Another sizable background component comes from 
isolated electrons, mainly originating from W and Z de- 
cays, which look similar to photons from the calorimetric 
point of view. The subtraction of such a contamination 
is addressed in Section IVIBI 

The background due to cosmic rays and to beam-gas 
collisions has been studied on dedicated data sets, se- 



lected by special triggers. Its impact is found to be neg- 
ligible. 



A. Jet background 

The jet background is due to photon-jet and di-jet final 
states. This section describes three methods, all based on 
the isolation transverse energy, E}fP , aiming to separate 
the TITI sample into four categories: 

AtTiti ^titi ^titi ^titi 

according to their physical final states — 7J and J7 dif- 
fer by the jet faking respectively the sub-leading or the 
leading photon candidate. The signal yield N'^F^'^ is eval- 
uated in bins of the three observables 771-77, -Pt,77, ^4>-i-yi 
as in Figure |3l Due to the dominant back-to-back topol- 
ogy of di-photon events, the kinematic selection produces 
a turn-on in the distribution of the di-photon invariant 
mass, at 77177 ~ 2i?™* (i?^"* = 16 GeV being the applied 
cut on the photon transverse energy), followed by the 
usual decrease typical of the continuum processes. The 
region at lower m^^ is populated by di-photon events 
with low A(/). 



The excess in the mass bin 80 < 






< 100 GeV, 



due to a contamination of electrons from Z-decays, is 
addressed in Section IVIBI 



From the evaluation of the background yields (A^^'-'^'- + 
iVj'F''^' and iVjl'''^'), the average fractions of photon-jet 
and di-jet events in the TITI sample are ^ 26% and 
~ 9% respectively. 

The three results shown in Figure |3] are compatible. 
This suggests that there are no hidden biases induced 
by the analyses. However, the three measures cannot be 
combined, as all make use of the same quantities — E}^° 
and shower shapes — and use the non-TiGHT background 
control region, so they may have correlations. None of 
the methods has striking advantages with respect to the 
others, and the systematic uncertainties are comparable. 
The "event weighting" method (IVI All) is used for the 
cross-section evaluation, since it provides event weights 
that are also useful in the event efficiency evaluation, and 
its sources of systematic uncertainties are independent of 
those related to the signal modelling and reconstruction. 



1. Event weighting 

Each event satisfying the tight selection on both pho- 
tons (sample TT) is classified according to whether the 
photons pass or fail the isolation requirement, resulting 
in a PP, PF, FP, or FF classification. These are trans- 
lated into four event weights VF77, W7J, iyj7, Wjj, which 
describe how likely the event is to belong to each of the 
four final states. A similar approach has already been 
used by the D0 |5| and CDF [6] collaborations. 
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FIG. 3. Differential 77 yields in the TITI sample {N^^'^^), 
as a function of the three observables m-yy, pT,~i-y, ^4>-tiy ob- 
tained with the three methods. In each bin, the yield is di- 
vided by the bin width. The vertical error bars display the 
total errors, accounting for both the statistical uncertainties 
and the systematic effects. The points are artificially shifted 
horizontally, to better display the three results. 



The connection between the pass/fail outcome and the 
weights, for the fc-th event, is: 



s 



S- 



pp 

(fe) 

PF 

(fc) 

FP 

(fc) 



= ^^^) 



\ 4'f^ / 



J7 



\ 



(1) 



If applied to a large number of events, the quantities Sxy 
would be the fractions of events satisfying each pass/fail 



classification, and the weights would be the fractions of 
events belonging to the four different final states. On an 
event- by- event approach, S\-y ^'^s boolean status vari- 
ables (e.g. for an event where both candidates are iso- 



lated, 4p == 1 and 4f 



4p = 4f = 0)- The 



quantity E'^*''' is a 4x4 matrix, whose coefficients give 
the probability that a given final state produces a cer- 
tain pass/fail status. If there were no correlation between 
the isolation transverse energies of the two candidates, it 
would have the form: 
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(2) 

where e^ and fi(i — l,2 for the leading/sub- leading can- 
didate) are the probabilities that a signal or a fake photon 
respectively pass the isolation cut. These are obtained 
from the E^° distributions extracted from collision data, 
as described in Section IV Al The value of e is essentially 
independent of E^ and changes with rj'^, ranging between 
80% and 95%. The value of / depends on both Et and 
ry and takes values between 20% and 40%. Given such 
dependence on the kinematics, the matrix E^*^' is also 
evaluated for each event. 

Due to the presence of correlation, the matrix coeffi- 
cients in equation ^ actually involve conditional prob- 
abilities, depending on the pass/fail status of the other 
candidate (tag) of the pair. For instance, the first two 
coefficients in the last column become: 



/i(l-/2)^ 



/f /2 + fifi 

/f(l-/2)+/l(l-/2^) 



where the superscripts P and F denote the pass/fail sta- 
tus of the tag. The ambiguity in the choice of the tag is 
solved by taking both choices and averaging them. All 

the conditional (e^ ' , /j ' ) probabilities are derived from 
collision data, as discussed in Section FV CI 

The signal yield in the TITI sample can be computed 
as a sum of weights running over all events in the TT 
sample: 



N, 



TITI 



= E ^^'^ ± 



,(k)] 



fceTT 



fcGTT 



where the weight w^''^ for the fc-th event is: 



(3) 



(4) 



and the sum over fc is carried out on the events in a given 
bin of the variable of interest {m^f^, Pt,77, ^4'-n)- The 
result is shown in Figure [3l by the solid circles. 

The main sources of systematic errors are: (?) the 
definition of the non-TiGHT control sample: ^ga- , (w) 
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FIG. 4. Projections of the 2-dimensional PDF fit on trans- 
verse isolation energies of the two photon candidates: leading 
photon (top) and sub- leading photon (bottom). Solid circles 
represent the observed data. The continuous curve is the fit 
result, while the dashed-dotted curve shows the 77 compo- 
nent. The dashed line represents the background component 
of the leading and sub-leading photon sample, respectively 



the normalization of the non-TiGHT sample: ^27' (*'*) 
the statistics used to compute the E}^" distributions, 
and hence the precision of the matrix coefficients: ±9%. 
Effects («) and (ii) are estimated as explained in Sec- 
tion |VAl Effect {in) is quantified by increasing and de- 
creasing the e, / parameters by their statistical errors, 
and recomputing the signal yield: the variations are then 
added in quadrature. 



2. Two-dimensional fit 

From all the di-photon events satisfying the tight 
selection (sample TT), the observed 2-dimensional dis- 
tribution F°^''{Ei^°^,E'r^°2) of the isolation energies of 
the leading and sub-leading photons is built. Then, 
a linear combination of four unbinned probability den- 
sity functions (PDFs), F^j, Fjj, Fjj, Fjj, describing the 2- 



dimensional distributions of the four final states, is fit to 
the observed distribution. For the 77, 7J, J7 final states, 
the correlation between £^x°i ^^^ -'^■¥'2 i^ assumed to be 
negligible, therefore the 2-dimensional PDFs are factor- 
ized into the leading and sub-leading PDFs. The leading 
and sub-leading photon PDFs F-y^ , F^^ are obtained from 
the electron extrapolation, as described in Section IV Bl 
The background PDF Fj^ for 7J events is obtained from 
the non-TlGHT sample on the sub-leading candidate, for 
events where the leading candidate satisfies the tight 
selection. The background PDF Fj^ for J7 events is ob- 
tained in a similar way. Both background PDFs are 
then smoothed with empirical parametric functions. The 
PDF for jj events cannot be factorized, due to the sizable 
correlation between the two candidates. Therefore, a 2- 
dimensional PDF is directly extracted from events where 
both candidates belong to the non-TiGHT sample, then 
smoothed. 

The four yields in the TT sample come from an ex- 
tended maximum likelihood fit of: 



N'-'-F°''%E'r^°^,E'r^°. 



7V,V^7l(£^T°l)^7.(^Tr2) 

iVj7Fj,(F'r-)F,,(i?^-2) 

■f^jj(,^T,l'-'^T,2J • 



JJ 



Figure 2] shows the fit result for the full TT data set. 

The yields in the TITI sample are evaluated by mul- 
tiplying N'F^ by the integral of the 2-dimensional sig- 
nal PDF in the region defined by El^\ < 3 GeV and 
i?^°2 < 3 GeV. The procedure is applied to the events be- 
longing to each bin of the observables ruj-y, Pt,77, A^^.^. 
The result is displayed in Figure[31 by the open triangles. 

The main sources of systematic uncertainties are: (i) 
definitionof the non-TiGHT control sample: Iq^ , (*J) sig- 
nal composition: ±8%; (iii) effect of material knowledge 
on signal: q^ , (iv) signal PDF parameters: ±0.7%; (v) 
jet PDF parameters: ±1.2%; (vi) di-jet PDF parameters: 
±1%; (vii) signal contamination in the non- tight sam- 
ple: Qw ■ Effect (i) is estimated by changing the num 



ber of released strips cuts, as explained in Section IV Al 
Effect (m) has been estimated by artificially setting the 
fraction of fragmentation photons to 0% or to 100%. Ef- 
fect {Hi) has been quantified by repeating the e — > 7 
extrapolation based on Monte Carlo samples with a dis- 
torted geometry. Effects {iv, v) have been estimated by 
randomly varying the parameters of the smoothing func- 
tions, within their covariance ellipsoid, and repeating the 
2-dimensional fit. Effect {vi) has been estimated by ran- 
domly extracting a set of {E}^°^, E}^°2) pairs, compara- 
ble to the experimental statistics, from the smoothed Fjj 
PDF, then re-smoothing the obtained distribution and 
repeating the 2-dimensional fit. Effect {vii) has been es- 
timated by taking the signal contamination from simula- 
tion — neglected when computing the central value. 
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FIG. 5. Schematic representation of the two-dimensional side- 
band method. The top plane displays the isolation (s-axis) 
and TIGHT identification (y-axis) criteria for the classification 
of the leading photon candidate. When the leading photon 
belongs to region A, the same classification is applied to the 
sub-leading photon, as described by the bottom plane. 



3. Isolation vs identification sideband counting 
(2D -sidebands) 

This method has been used in ATLAS in the inclusive 
photon cross-section measurement [lOi] and in the back- 
ground decomposition in the search for the Higgs boson 
decaying into two photons |19j . 

The base di-photon sample must fulfil the selection 
with the strips cuts released, defined by the union of 
TIGHT and non-TiGHT samples and here referred to as 
loose' (L'). The leading photons in the L'L' sample 
are divided into four categories A^ B, C, D, depending 
on whether they satisfy the tight selection and/or the 
isolation requirement — see Figure [S] (top). The signal 
region, defined by tight and isolated photons (TI), con- 
tains Na candidates, whereas the three control regions 
contain Nb, Nq, Nd candidates. Under the hypothesis 
that regions B, C, D are largely dominated by back- 
ground, and that the isolation energy of the background 



has little dependence on the tight selection (as discussed 
in Section fV Ap . the number of genuine leading photons 
N^^ in region A, coming from 77 and 7J final states, can 
be computed [l3| by solving the equation: 



NT- 



Na 



{Nb - ci7V7) 



Nc - C2NT 
Nd - C1C2N' 



^bkg 



(5) 

Here, ci and C2 are the signal fractions failing respec- 
tively the isolation requirement and the tight selection. 
The former is computed from the isolation distributions, 
as extracted in Section fV A[ the latter is evaluated from 
Monte Carlo simulation, after applying the corrections 
to adapt it to the experimental shower shapes distribu- 

Arbkg iirbkg 

tions [l3|. The parameter R^^'^ — bkg -.bicg measures 

the degree of correlation between the isolation energy 
and the photon selection in the background: it is set to 1 
to compute the central values, then varied according to 
the "di-jet-like" Monte Carlo prediction for systematic 
studies. 

When the leading candidate is in the TI region, 
the sub-leading one is tested, and four categories 
A', B', C", D' are defined, as in the case of the lead- 
ing candidate — see Figure [5] (bottom). The number of 
genuine sub-leading photons N'^^^ , due to 77 and J7 fi- 
nal states, is computed by solving an equation analogous 
to ®. 

N^ and N'j^^^^ are related to the yields by: 



NT- 



N 



N'/'^ ^ N 



TITI 

77 

e' 

TITI 



^TITI 

7J 



77 



N 



f 

TITI 

J7 



d+c' )(i+c' ) ^^ ^^*^ probability that a sub- 
leading photon satisfies the tight selection and isolation 
requirement, while /' is the analogous probability for a 
jet faking a sub-leading photon. The di-photon yield is 
therefore computed as: 



N 



TITI 

77 



e' («/' 



NT + {a- 



m'/'] 



{a~l)e' + af' 



(6) 



and /' can be computed from the observed quantities to 



be/' = 



N'.-N' 



The parameter a is defined as the 



fraction of photon-jet events in which the jet fakes the 
leading photon, a = ^titiImtiti ; whose value is taken 
from the Pythia photon-jet simulation. 



The counts Na, Nb, Nc, Nd, N'^, N'g, N^, N'^, 



and 



hence the yield, can be computed for all events entering 
a given bin of m^-y, pt,77, A^-y.^. The result is displayed 
in Figure [21 by the open squares. 

The main source of systematic error is the definition 
of the non-TiGHT sample: it induces an error of I^qV ■ 
The other effects come from the uncertainties of the pa- 
rameters entering equation ^. The main effects come 
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FIG. 6. Electron background subtraction as a function of 
m^.y. The top plot displays the impurity, overall and for 
the 76 and ee separately. The bottom plot shows the di- 
photon yield before (open squares) and after (solid circles) the 
electron background subtraction. The points are artificially 
shifted horizontally, to better display the different values. 



electrons satisfying a tight selection and the calorimet- 
ric isolation E^° < 3 GeV are considered, and electrons 
are counted only if they are not reconstructed at the same 
time as photons. Such counts are related to the actual 
underlying yields iV*™°, A^.^™'', iV*™°, defined as the num- 
ber of reconstructed final states where both particles are 
correctly classified. Introducing the ratio /e->7 = jf-^^ 
between genuine electrons that are wrongly and correctly 
classified, and likewise f-y^e — Af^' ^oi genuine photons, 
the relationship between the N and 7V*''"° quantities is 
described by the following linear system: 



1 /e^T (/e^- 






which can be solved for the unknown 7V'™°. 



(7) 



77 



The value of fe~^j is extracted from collision data, as 

fe^-y — 2N° ' from events with an invariant mass within 
±5 GeV of the Z mass. The continuum background 
is removed using symmetric sidebands. The result is 
/e^7 = 0.112 ± 0.005(stat) ± 0.003(syst), where the sys- 
tematic error comes from variations of the mass window 
and of the sidebands. This method has been tested on 
"di-jet-like" and Z -> e+e^ Monte Carlo samples and 
shown to be unbiased. The value of f^^e is taken from 
the "di-jet-like" Monte Carlo: /^^e = 0.0077. To ac- 
count for imperfect modelling, this value has also been 
set to 0, or to three times the nominal value, and the re- 
sulting variations are considered as a source of systematic 
error. 

The electron contamination is estimated for each bin 
of TO.^^, pt,77 and A0^.^, and subtracted from the di- 
photon yield. The result, as a function of m^^, is shown 
in Figure |B1 The fractional contamination as a function 
of pt,77 E^nd A0-y.^ is rather flat, amounting to ^ 5%. 



from: (?) variation of c'j^ : ±4%; [it) variation of a: ±3%; 
(m) variations of i^'^'^s^ i^'^kg. +o%_^^^ The variations of 

ci, C2, c'2 have negligible impact. 



B. Electron background 

Background from isolated electrons contaminates 
mostly the selected converted photon sample. The con- 
tamination in the di-photon analysis comes from several 
physical channels: (i) e+e" final states from Drell-Yan 
processes, Z — >• e+e" decay, W^W~ — )■ e'^e~vv] {ii) 
76^ final states from di-boson production, e.g. ^W^ — >■ 
76^1^, 7Z — >■ 7e'*"e~. The effect of the Z — >■ e+e" 
contamination is visible in Figure [3] in the mass bin 



< TO77 < 100 GeV. 



Rather than quantifying each physical process sepa- 
rately, a global approach is chosen. The events recon- 
structed with 77, 7e and ee final states are counted, thus 
obtaining counts ^^^7, iV^e and A^ee- Only photons and 



VII. EFFICIENCIES AND UNFOLDING 

The signal is defined as a di-photon final state, which 
must satisfy precise kinematic cuts (referred to as "fidu- 
cial acceptance" ) : 

• both photons must have a transverse momentum 
p^ > 16 GeV and must be in the pseudorapidity 
acceptance \Tf\ < 2.37, with the exclusion of the 
region 1.37 < \7f\ < 1.52; 

• the separation between the two photons must be 
Ai?77 = VivJ - VIY + W - 4>IY > 0.4; 

• both photons must be isolated, i.e. the transverse 
energy fiow i5^° due to interacting particles 
in a cone of angular radius R < 0.4 must be 



£'■ 



3(part) 



< 4 GeV. 



These kinematic cuts define a phase space similar to the 
experimental selection described in Section Hill In partic- 
ular, the requirement on El^° ^^^ has been introduced 
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to match approximately the experimental cut on E}^°. E^ > 16 GeV. Finally, di-photon TITI events with 

the sub-leading photon selected by a high-level photon 
trigger are used to compute the efficiency of the di- 
photon 15 GeV-threshold high-level trigger, obtaining 
£2 — 99A^1'q%. The overall efficiency of the trigger is 
therefore e*"59'='~ = £06162 = (98.5t?;|5 ± 1.0)%. The first 
uncertainty is statistical, the second is systematic and 
accounts for the contamination of photon-jet and di-jet 
events in the selected sample. 



The value of £'^° is corrected for the ambient en- 

ergy, similarly to what is done for E}^°. From studies on 
a Pythia di-photon Monte Carlo sample, there is a high 
correlation between the two variables, and E^° = 3 GeV 

corresponds to E^r^°^'^^'^ ' ~ 4 GeV. 

A significant number of di-photon events lying outside 
the fiducial acceptance pass the experimental selection 
because of resolution effects: these are referred to as "be- 
low threshold" (BT) events. 

The background subtraction provides the di-photon 
signal yields for events passing all selections (TITI). 
Such yields are called N^^"^^, where the index i flags 
the bins of the reconstructed observable X'^°'' under con- 
sideration (X being ruj^, PT,-yy^ ^'P-y-y)- The relationship 
between TV^'''^' and the true yields Ua {a being the bin 
index of the true value X'''"'') is: 



m^ii^L 



^TITI 

BT\ 



^trigger ^TT j^ll 



= 2J Mjac„ 






(8) 
(9) 



where iV" is the number of reconstructed isolated di- 
photon events in the i-th bin, and: 

• e^"99er jg ^jj^g trigger efficiency, computed for events 
where both photons satisfy the tight identification 
and the calorimetric isolation; 

• e^"^ is the efficiency of the tight identification, for 
events where both photons satisfy the calorimetric 
isolation; 

• ff^ is the fraction of "below-threshold" events; 

• Mia is a "migration probability" , i.e. the proba- 
bility that an event with x'"^"" in bin-a is recon- 
structed with X™^ in bin-i; 

• e^^ accounts for both the reconstruction efficiency 
and the acceptance of the experimental cuts (kine- 
matics and calorimetric isolation). 



A. Trigger efficiency 

The trigger efficiency is computed from collision data, 
for events containing two reconstructed photons with 
transverse energy E^ > 16 GeV, both satisfying the 
tight identification and the calorimetric isolation re- 
quirement (TITI). The computation is done in three 
steps. 

First, a level- 1 e/7 trigger with an energy threshold 
of 5 GeV is studied: its efficiency, for reconstructed TI 
photons, is measured on an inclusive set of minimum- 
bias events: for E^^ > 16 GeV it is eo = lOCOtH^ 
— therefore such a trigger does not bias the sample. 
Next, a high-level photon trigger with a 15 GeV thresh- 
old is studied, for reconstructed TI photons selected by 
the level-1 trigger: its efficiency is ei = 99.ltQ;4% for 



B. Identification efficiency 



The photon tight identification efficiency e ' , for 
photon candidates satisfying the isolation cut E^° < 
3 GeV, is computed as described in Ref [Ifll , as a function 
of 77''' and E^. The efficiency is determined by applying 
the TIGHT selection to a Monte Carlo photon sample, 
where the shower shape variables have been shifted to 
better reproduce the observed distributions. The shift 
factors are obtained by comparing the shower shapes 
of photon candidates from a "di-jet-like" Monte Carlo 
sample to those observed in collision data. To enhance 
the photon component in the sample — otherwise over- 
whelmed by the jet background — only the photon candi- 
dates satisfying the tight selection are considered. This 
procedure does not bias the bulk of the distribution un- 
der test appreciably, since the cuts have been tuned to 
reject only the tails of the photons' distributions. How- 
ever, to check the systematic effect due to the selection, 
the shift factors are also recomputed applying the LOOSE 
selection. 

Compared to Ref [10], the photon identification cuts 
have been re-optimized to reduce the systematic er- 
rors, and converted and unconverted photons treated 
separately. The photon identification efficiency is rj'^- 
dependent, and increases with E^, ranging from ~ 60% 
for 16 < e:^ < 20 GeV, to > 90% for £;^ > 100 GeV. 
The overall systematic error is between 2% and 10%, the 
higher values being applicable at lower E^ and for con- 
verted photons. The main sources of systematic uncer- 
tainty are: (i) the systematic error on the shift factors; 
(ii) the knowledge of the detector material; {in) the fail- 
ure to detect a conversion, therefore applying the wrong 
TIGHT identification. 

Rather than computing an event- level identification ef- 
ficiency for each bin of each observable, the photon ef- 
ficiency can be naturally accommodated into the event 
weights described in Section IVIAll by dividing the 
weight w'-'^^ of equation (JH) by the product of the two 
photon efficiencies: 



N. 



II 



(k) 



V 



(fc) 



(10) 



where the sum is extended over all events in the TT sam- 
ple and in the i-th bin. Here the identification efficien- 
cies of the two photons are assumed to be uncorrelated 
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— which is ensured by the separation cut AR > 0.4, and 
by the binning in rj'^ and E^. 

The event efficiency, ef^ — —krn-, is essentially flat 

at '^ 60% in A0^^, and increases with m^^ and pt,77, 
ranging from ~ 55% to ~ 80%. Its total systematic error 
is '^ 10%, rather uniform over the m^j, PT,'y-y, A0-y^ 
ranges. 



C. Reconstruction, acceptance, isolation and 
unfolding 

The efficiency e^^ accounts for both the reconstruc- 
tion efficiency and the acceptance of the experimental 
selection. It is computed for each bin of X'''"°, with 
Monte Carlo di-photon events generated with Pythia in 
the fiducial acceptance, as the fraction of events where 
both photons are reconstructed, pass the acceptance cuts 
and the calorimetric isolation. The value of e^"^ ranges 
between 50% and 60%. The two main sources of ineffi- 
ciency are the local ECAL read-out failures (^ —18%) 
and the calorimetric isolation (^ —20%). 

The energy scale differences between Monte Carlo and 
collision data — calibrated on Z — >■ e~^e~ events — are 
taken into account. The uncertainties on the energy 
scale and resolution are propagated as systematic errors 
through the evaluation: the former gives an effect be- 
tween +3% and —1% on the signal rate, while the latter 
has negligible impact. 

In Monte Carlo, the calorimetric isolation energy, E^^", 
needs to be corrected to match that observed in colli- 
sion data. The correction is optimized on tight pho- 
tons, for which the background contamination can be 
removed (see Section fVAj) . then it is applied to all pho- 
tons in the Monte Carlo sample. The E!^° difference 
observed between Monte Carlo simulation and collision 
data may be entirely due to inaccurate GEANT4/detector 
modelling, or it can also be a consequence of the physical 
model in the generator (e.g. kinematics, fragmentation, 
hadronization) . From the comparison between collision 
data and simulation, the two effects cannot be disentan- 
gled. To compute the central values of the results, the dif- 
ference between simulation and collision data is assumed 
to be entirely due to the detector simulation. As a cross- 
check, the opposite case is assumed: that the difference 
is entirely due to the generator model. In this case, the 
particle-level isolation El^^^^'^ ' should also be corrected, 

using the i5^° -^ E}f° relationship described by the 

detector simulation. This modifies the definition of fidu- 
cial acceptance, and hence the values of e^^, resulting in 
a cross-section variation of ~ —7%, which is handled as 
an asymmetric systematic uncertainty. 

The fraction of events "below threshold" , fP"^, is com- 
puted from the same Pythia signal Monte Carlo sample, 
for each bin of X^°^. Its value is maximum [^ 12%) for 
m-y-y about twice the -Et cut, and decreases to values 
< 5% for m^^ > 50 GeV. 



The "migration matrix", Mia, is filled with Pythia 
Monte Carlo di-photon events in the fiducial acceptance, 
that are reconstructed, pass the acceptance cuts and the 
calorimetric isolation. The inversion of this matrix is per- 
formed with an unfolding technique, based on Bayesian 
iterations [20| . The systematic uncertainties of the pro- 
cedure have been estimated with a large number of toy 
datasets and found to be negligible. The result has also 
been tested to be independent of the initial ( "prior" ) dis- 
tributions. Moreover, it has been checked that a simpler 
bin-by-bin unfolding yields compatible results. 

As the evaluation of e^"^, f^^, Mia may strongly de- 
pend on the simulation modelling, two additional Monte 
Carlo samples have been used, the first with more mate- 
rial modelled in front of the calorimeter, and the second 
with a different generator (Sherpa): the differences on 
the computed signal rates are ~ -1-10% and < +5% re- 
spectively, and are treated as systematic errors. 



VIII. CROSS-SECTION MEASUREMENT 

The di-photon production cross-section is evaluated 
from the corrected binned yields Ua , divided by the inte- 
grated luminosity / Mt = (37.2 ± 1.3) pb'^ Q. The re- 
sults are presented as differential cross-sections, as func- 
tions of the three observables m-^-y, pt,77, A^^.^, for a 
phase space defined by the fiducial acceptance cuts in 
Section I VIII In Table HI the differential cross-section is 
quoted for each bin, with its statistical and systematic 
uncertainty. In Table [III all the considered sources of 
systematic errors are listed separately. 

The experimental measurement is compared with 
theoretical predictions from the DIPHOX 21] and 
ResBos '22^ NLO generators in Figures H 1] and HI 
The DIPHOX and ResBos evaluation has been carried 
out using the NLO fragmentation function |2J a,nd the 
CTEQ6.6 parton density function (PDF) setj^- The 
fragmentation, normalization and factorization scales are 
set equal to m^-y. The same fiducial acceptance cuts in- 
troduced in the signal definition (Section IVII[) are ap- 
plied. Since neither generator models the hadronization, 
it is not possible to apply a requirement on i5^° : 

the closest isolation variable available in such generators 
is the "partonic isolation", therefore this is required to 
be less then 4 GeV. The computed cross-section shows 
a weak dependence on the partonic isolation cut: mov- 
ing it to 2 GcV or 6 GcV produces variations within 5%, 
smaller than the theoretical systematic errors. 

The theory uncertainty error bands come from scale 
and PDF uncertainties evaluated from DIPHOX: (i) vari- 
ation of rcnormalization, fragmentation and factorization 
scales: each is varied to ^m.77 and 2m.y^, and the en- 
velope of all variations is assumed as a systematic er- 
ror; (ii) variation of the eigenvalues of the PDFs: each 
is varied by ilcr, and positive/negative variations are 
summed in quadrature separately. As an alternative, the 
MSTW 2008 PDF set has been used: the difference with 
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FIG. 7. Differential cross-section Aa jAm-^r^ of di-photon pro- 
duction. The solid circles display the experimental values, the 
hatched bands display the NLO computations by DIPHOX 
and ResBos. The bottom panels show the relative difference 
between the measurements and the NLO predictions. The 
data/theory point in the bin < m^f-^ < 30 GeV lies above 
the frames. 



FIG. 8. Differential cross-section da/dpT,-yy of di-photon 
production. The solid circles display the experimental val- 
ues, the hatched bands display the NLO computations by 
DIPHOX and ResBos. The bottom panels show the relative 
difference between the measurements and the NLO predic- 
tions. The data point in the bin 150 < Pt,77 < 200 GeV in 
the main panel lies below the frame. 



IX. CONCLUSIONS 



respect to CTEQ6.6 is an overall increase by ~ 10%, 
which is covered by the CTEQ6.6 total systematic error. 



The measured distribution of da / dA(j)j.y (Figure [5]) is 
clearly broader than the DIPHOX and ResBos predic- 
tions: more photon pairs are seen in data at low Acj}^^ 
values, while the theoretical predictions favour a larger 
back-to-back production (Ac/i^-y ~ tt). This result is qual- 
itatively in agreement with previous measurements at the 
Tevatron [^, Q . The distribution of da/dm^-y (Figure [7]) 
agrees within the assigned uncertainties with both the 
DIPHOX and ResBos predictions, apart from the region 
rrij^ < 2£;^"* {E^''^ = 16 GeV being the applied cut on 
the photon transverse momenta): as this region is popu- 
lated by events with small A(j)j-y , the poor quality of the 
predictions can be related to the discrepancy observed in 
the A(j)jj distribution. The result for da/dpr.-yj (Fig- 
ure [5]) is in agreement with both DIPHOX and ResBos: 
the maximum deviation, about 2a, is observed in the re- 
gion 50 < pt,77 < 60 GeV. 



This paper describes the measurement of the pro- 
duction cross-section of isolated di-photon final states 
in proton-proton collisions, at a centre-of-mass energy 
y/s = 7 TeV, with the ATLAS experiment. The full data 
sample collected in 2010, corresponding to an integrated 
luminosity of 37.2 ± 1.3 pb~^, has been analysed. 

The selected sample consists of 2022 candidate events 
containing two reconstructed photons, with transverse 
momenta px > 16 GeV and satisfying tight identifi- 
cation and isolation requirements. All the background 
sources have been investigated with data-driven tech- 
niques and subtracted. The main background source, 
due to hadronic jets in photon-jet and di-jet events, has 
been estimated with three computationally independent 
analyses, all based on shower shape variables and isola- 
tion, which give compatible results. The background due 
to isolated electrons from W and Z decays is estimated 
with collision data, from the proportions of observed ee, 
76 and 77 final states, in the Z-mass region and else- 
where. 

The result is presented in terms of differential cross- 
sections as functions of three observables: the invariant 
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FIG. 9. Differential cross-section da/dA^T,^ of di-photon 
production. The solid circles display the experimental val- 
ues, the hatched bands display the NLO computations by 
DIPHOX and ResBos. The bottom panels show the relative 
difference between the measurements and the NLO predic- 
tions. 



mass m^^^ the total transverse momentuni pt,77, and 
the azimuthal separation ts.<^^^ of the photon pair. The 
experimental results are compared with NLO predictions 
obtained with DIPHOX and ResBos generators. The ob- 
served spectrum of da/dA^.^-^ is broader than the NLO 
predictions. The distribution of da/dm^-y is in good 
agreement with both the DIPHOX and ResBos predic- 
tions, apart from the low mass region. The result for 
dcr/dpT,77 is generally well described by DIPHOX and 
ResBos. 



TABLE I. Binned differential cross-sections da/dm-y-,, 
da/dpT,77, dcr/dAt/)^.^ for di-photon production. For each 
bin, the differential cross-section is quoted with its statisti- 
cal and systematic uncertainties (symmetric and asymmetric, 
respectively). Values quoted as 0.000 are actually less than 
0.0005 in absolute value. 



[GeV] dajdm.,^ [pb/GeV] 





30 
40 
50 
60 
70 



30 
40 
50 
60 
70 
80 



80 - 100 
100 - 150 
150 - 200 



0.20 ± 0.05 


+0.05 
-0.03 


1.8 ± 0.3 


+0.4 
-0.3 


2.3 ± 0.3 


+0.6 
-0.4 


1.83 ± 0.24 


+0.36 
-0.28 


0.74 ± 0.17 


+0.19 
-0.13 


0.45 ± 0.15 


+0.11 
-0.09 


0.40 ± 0.06 


+0.08 
-0.08 


0.079 ± 0.022 


+0.025 
-0.025 



0.026 ± 0.009 to.Sw 



Pt,77 [GeV] d(T/dpT,77 [pb/GeV] 



+0.9 

-0.6 

+0.5 

-0.4 

+0.28 

-0.24 

+0.21 

-0.14 

+0.10 

-0.09 

+0.09 

-0.05 

+0.03 

-0.03 

+0.009 

-0.010 

+0.003 

-0.002 

+0.000 

-0.000 



A.(j).,^ [rad] da/dA(l>~/-, [pb/rad] 



+ 1.5 

-1.1 

+2.5 
-1.9 
+6 
-4 

+ 12 

-9 

+24 

-17 

+36 

-29 






- 10 


4.5 


± 0.4 


10 


- 20 


2.2 


±0.3 


20 


- 30 


0.94 


± 0.22 


30 


- 40 


0.62 


± 0.16 


40 


- 50 


0.26 


± 0.10 


50 


- 60 


0.36 


± 0.09 


60 


- 80 


0.06 


±0.03 


80 


- 100 


0.048 


± 0.019 


100 


- 150 


0.003 


± 0.004 


150 


- 200 


0.000 


± 0.002 



0.00 - 


1.00 


4.9 ± 1.1 


1.00 - 


2.00 


8.9 ± 1.8 


2.00 - 


2.50 


24 ± 4 


2.50 - 


2.80 


56 ± 8 


2.80 - 


3.00 


121 ± 13 


3.00 ~ 


3.14 


173 ± 16 
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TABLE II. Breakdown of the total cross-section systematic uncertainty, for each bin of m~f^, PT,-fy and Acfi^-,. The meaning 
of each column is as follows: "T" is the definition of the non-TlGHT control sample; "I" is the choice of the E^° region used 
to normalize the non-TlGHT sample; "matrix" refers to the statistical uncertainty of the matrix coefficients used by the event 
weighting; "e — > 7" is the total systematic coming from the electron fake rate; "ID" is the overall uncertainty coming from the 
method used to derive the identification efficiency; "material" is the effect of introducing a detector description with distorted 
material distribution; "generator" shows the variation due to the usage of a difi'erent generator (Sherpa instead of Pythia); 
"(Tb" and "i5-scale" are due to uncertainties on energy resolution and scale; "^^^"(p''''')" jg the effect of smearing the particle- 



liso(part) 

'T ; 

than 0.0005 in absolute value. 



level isolation i?.^°^'""^'' ; "J Ldt" is the effect due to the total luminosity uncertainty. Values quoted as 0.000 are actually less 



rrij^ [GeV] 



matrix e 



ID material generator ge E-scale E\ 



so(part) 



j-LAt 



0-30 
30 - 40 
40 - 50 
50 - 60 
60 - 70 
70 - 80 
80 - 100 
100 - 150 
150 - 200 



+0.03 
-0.01 


+0.000 
-0.005 


+0.021 
-0.022 


+0.002 
-0.002 


+0.020 
-0.017 


+0.021 
-0.000 


+0.03 
-0.00 


+0.001 
-0.000 


+0.006 
-0.002 


+0.000 
-0.010 


+0.007 
-0.007 


+0.17 
-0.09 


+0.00 
-0.05 


+0.13 
-0.13 


+0.008 
-0.008 


+0.22 
-0.18 


+0.3 
-0.0 


+0.014 
-0.000 


+0.003 
-0.000 


+0.04 
-0.03 


+0.00 
-0.09 


+0.06 
-0.06 


+0.3 
-0.1 


+0.00 
-0.06 


+0.19 
-0.19 


+0.008 
-0.008 


+0.24 
-0.20 


+0.3 
-0.0 


+0.11 
-0.00 


+0.024 
-0.000 


+0.09 
-0.03 


+0.00 
-0.17 


+0.08 
-0.08 


+0.20 
-0.13 


+0.00 
-0.04 


+0.14 
-0.14 


+0.007 
-0.007 


+0.15 
-0.13 


+0.19 
-0.00 


+0.05 
-0.00 


+0.003 
-0.000 


+0.06 
-0.03 


+0.00 
-0.13 


+0.06 
-0.06 


+0.14 
-0.06 


+0.001 
-0.016 


+0.09 
-0.09 


+0.004 
-0.004 


+0.05 
-0.04 


+0.07 
-0.00 


+0.04 
-0.00 


+0.007 
-0.000 


+0.03 
-0.02 


+0.00 
-0.05 


+0.03 
-0.03 


+0.06 
-0.06 


+ 0.000 
-0.007 


+ 0.05 
-0.06 


+0.003 
-0.003 


+0.03 
-0.03 


+0.07 
-0.00 


+0.03 
-0.00 


+0.002 
-0.001 


+0.009 
-0.002 


+0.00 
-0.03 


+0.015 
-0.015 


+0.04 
-0.05 


+0.000 
-0.005 


+0.04 
-0.04 


+0.019 
-0.019 


+0.03 
-0.02 


+0.04 
-0.00 


+0.012 
-0.000 


+0.004 
-0.000 


+0.013 
-0.003 


+0.00 
-0.03 


+0.013 
-0.013 


+0.019 
-0.016 


+0.001 
-0.001 


+0.015 
-0.018 


+0.001 
-0.001 


+0.004 
-0.003 


+0.002 
-0.001 


+0.004 
-0.000 


+0.000 
-0.001 


+0.002 
-0.003 


+0.000 
-0.005 


+0.003 
-0.003 


+0.002 
-0.002 


+0.000 
-0.000 


+0.003 
-0.003 


+0.000 
-0.000 


+0.002 
-0.001 


+0.004 
-0.000 


+0.001 
-0.000 


+0.000 
-0.000 


+0.001 
-0.000 


+0.000 
-0.002 


+0.001 
-0.001 



Pt,77 [GeV] 



matrix e - 



ID material generator oe B-scalc E\ 



so(part) 



/Ldt 



0-10 
10 - 20 
20 - 30 
30 - 40 
40 - 50 
50 - 60 
60 - 80 
80 - 100 
100 - 150 
150 - 200 



+0.3 

-0.2 


+0.00 
-0.09 


+0.3 
-0.3 


+0.03 
-0.03 


+0.4 
-0.4 


+0.6 
-0.0 


+0.10 
-0.00 


+0.03 
-0.00 


+0.12 
-0.05 


+0.0 
-0.3 


+0.15 
-0.15 


+0.3 
-0.2 


+0.00 
-0.05 


+0.21 
-0.22 


+0.015 
-0.015 


+0.20 
-0.17 


+0.21 
-0.00 


+0.11 
-0.00 


+0.001 
-0.001 


+0.06 
-0.03 


+0.00 
-0.15 


+0.08 
-0.08 


+0.21 
-0.16 


+0.000 
-0.025 


+0.13 
-0.14 


+0.008 
-0.008 


+0.07 
-0.06 


+0.10 
-0.00 


+0.022 
-0.000 


+0.010 
-0.000 


+0.03 
-0.02 


+0.00 
-0.08 


+0.03 
-0.03 


+0.13 
-0.08 


+0.000 
-0.012 


+0.09 
-0.10 


+0.006 
-0.006 


+0.06 
-0.05 


+0.11 
-0.00 


+0.08 
-0.00 


+0.007 
-0.000 


+0.015 
-0.009 


+0.00 
-0.03 


+0.021 
-0.021 


+0.08 
-0.06 


+ 0.000 
-0.007 


+0.05 
-0.06 


+0.004 
-0.004 


+0.018 
-0.017 


+0.03 
-0.00 


+0.005 
-0.000 


+0.000 
-0.012 


+0.00 
-0.03 


+0.000 
-0.015 


+0.009 
-0.009 


+0.03 
-0.03 


+0.000 
-0.007 


+0.02 
-0.03 


+0.006 
-0.006 


+0.03 
-0.02 


+0.04 
-0.00 


+0.04 
-0.00 


+0.013 
-0.000 


+0.05 
-0.01 


+0.000 
-0.023 


+0.012 
-0.012 


+0.021 
-0.023 


+0.000 
-0.001 


+0.014 
-0.016 


+0.001 
-0.001 


+0.003 
-0.003 


+0.005 
-0.000 


+0.000 
-0.004 


+0.000 
-0.001 


+0.000 
-0.002 


+0.000 
-0.004 


+0.002 
-0.002 


+0.006 
-0.000 


+0.000 
-0.001 


+0.005 
-0.005 


+0.002 
-0.002 


+0.003 
-0.002 


+0.002 
-0.006 


+0.000 
-0.005 


+0.001 
-0.000 


+0.004 
-0.001 


+0.000 
-0.004 


+0.002 
-0.002 


+0.002 
-0.001 


+0.000 
-0.000 


+0.001 
-0.002 


+0.000 
-0.000 


+0.000 
-0.000 


+0.000 
-0.000 


+0.000 
-0.001 


+0.000 
-0.000 


+0.000 
-0.000 


+0.000 
-0.000 


+0.000 
-0.000 


+0.000 
-0.000 


+ 0.000 
-0.000 


+0.000 
-0.000 


+0.000 
-0.000 


+0.000 
-0.000 


+0.000 
-0.000 


+0.000 
-0.000 


+0.000 
-0.000 


+0.000 
-0.000 


+0.000 
-0.000 


+ 0.000 
-0.000 



A<^7 



ad] 



i matrix e — > 7 ID material generator oe 



E-scalc £^=°(P-') 



JLdt 



0.00 - 1.00 

1.00 - 2.00 

2.00 - 2.50 

2.50 - 2.80 

2.80 - 3.00 

3.00 - 3.14 



1.1 
0.5 


+0.00 
-0.14 


+ 0.8 
-0.8 


+0.05 
-0.05 


+0.4 
-0.4 


+0.4 
-0.0 


+0.3 
-0.0 


+0.000 
-0.017 


+0.14 
-0.08 


+0.0 
-0.3 


+0.17 
-0.17 


1.6 
1.0 


+ 0.0 
-0.3 


+ 1.2 
-1.2 


+0.07 
-0.07 


+0.8 
-0.7 


+ 1.0 
-0.0 


+0.5 
-0.0 


+0.023 
-0.000 


+0.23 
-0.10 


+0.0 
-0.5 


+0.3 
-0.3 


3 

2 


+0.0 
-0.4 


+2.2 
-2.3 


+0.17 
-0.17 


+2.2 
-1.8 


+3 
-0 


+ 1.5 
-0.0 


+0.10 
-0.00 


+0.6 
-0.4 


+0.0 
-1.3 


+0.8 
-0.8 


6 

5 


+0.0 
-1.3 


+5 
-5 


+0.4 
-0.4 


+5 
-4 


+6 
-0 


+0.3 
-0.0 


+0.4 
-0.0 


+ 1.8 
-1.0 


+0 
-4 


+ 1.9 
-1.9 


11 


+0 
-3 


+9 
-10 


+0.9 
-0.9 


+ 11 
-9 


+ 14 
-0 


+2.3 
-0.0 


+0.7 
-0.0 


+4 
-1 


+0 
-9 


+4 
-4 


19 
16 


+0 
-3 


+ 14 
-15 


+ 1.5 
-1.5 


+ 16 
-13 


+ 18 
-0 


+9 
-0 


+0.6 
-0.0 


+4 
-2 


+0 
-12 


+6 
-6 
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